Unclassified 


AD-A269  874 


REPORT  DOCUMENTS 


la.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

Za.  SECURITY  CLASSIFICATION  AUTHORITY 

2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

4.  PERFORM.NG  ORGANIZATION  REPORT  NUM8ER($) 

6a.  NAME  OF  PERFORMING  ORGANIZATION 

Calif.  Institute  of  Technology 

6b  OFFICE  SYMBOL 
(If  applicable) 

6c  ADDRESS  (Cfty,  Sfata,  and  ZIP  Code) 

Mail  Code  104-44 

Pasadena,  CA  91125 

8a.  NAME  OF  FUNDING /SPONSORING 

9b.  OFFICE  SYMBOL 

ORGANIZATION 

Office  of  Naval  Research 

(If  applicable) 

8c  ADDRESS  fCJfy.  State,  and  ZIP  Code) 

800  N.  Quincy  St. 

Arlington,  VA  22217-5000 

lb  RESTRfCTIVE  MARKINGS 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Office  of  Naval  Research  (Dr.  Dan  Harris) 


7b.  ADDRESS  (Ofy,  Sfa»a,  and  ZiPCodt) 

Code  3854 

Naval  Air  Warfare  Center 
China  Lake,  C 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


1 1 .  TITLE  (/nc/ude  Security  Clauification) 

Scaling  Laws  for  Diamond  Chemical  Vapor  Deposition.  I:  Diamond  Surface  Chemistry 


1Z.  PERSONAL  AUTHOR(S) 

David  G.  Goodwin 


13a.  TYPE  OF  REPORT 

technical 


16.  SUPPLEMENTARY  NOTATION 


13fa  TIME  COVERED 
FROM  TO 


14,  DATE  OF  REPORT  (Year,  Month,  Day)  15  PAGE  COUNT 

93/9/15  19 


COSATI  COOES 


GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverie  If  neceaary  and  identify  by  block  number) 


19.  ABSTRACT  {Continue  on  revene  If  necessary  and  Identify  by  block  number) 


A  simplified  model  of  the  gas-surface  chemistry  occurring  during  chamical  vapor 
deposition  of  diamond  thin  films  is  presented.  The  model  results  in  simple  scaling 
relations,  useful  for  process  scale-up  and  optimization,  for  growth  rate  and  defect 
density  in  terms  of  the  local  chemical  environment  at  the  substrate.  A  simple 
two-parameter  expression  for  growth  rate  is  obtained,  which  with  suitable  parameter 
choices  reproduces  the  results  of  more  detailed  mechanisms  and  experiment  over  two 
orders  of  magnitude  in  growth  rate.  The  defect  formation  model  suggests  that  the 
achievable  growth  rate  at  specified  defect  density  scales  approximately  quadratir-oTi 
with  the  atomic  hydrogen  concentration  at  the  substrate. 


C 


0  0 


93-2;^ 


I 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
GJunCLASSIFIEO/UNLIMITED  □  SAME  AS  RPT.  Q  OTIC  USERS 

21.  ABSTRACT  SEC 

_ Unclassi 

\  hi 

_ 

22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

David  G.  Goodwin 

22b  TELEPHONE  (Include  Area  Code) 

818-395-4249 

22c  OFFICE  SYMBOL 

Form  1473,  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


to  appear  in  Journal  of  Applied  Physics,  1  Dec  1993 


Scaling  Laws  for  Diamond  Chemical  Vapor  Deposition. 
I:  Diamond  Surface  Chemistry 


Accebio,.'  For 

NTI$  CRA&l 
UIIC  lAB 


D.  G.  Goodwin 

Division  of  Engineering  and  Applied  Science 
California  Institute  of  Technology 
Pasadena,  CA  91125 


t; 


By. 

Diht-ib 

Aionl  1 

Avaiidbiiily  Coaes  | 

Dist 

Avail  ant';  or 

Special 

A  simplified  model  of  the  gas-surface  chemistry  occurring  during  chemical  vapor 
deposition  of  diamond  thin  films  is  presented.  The  model  results  in  simple  scaling 
relations,  useful  for  process  scale-up  and  optimization,  for  growth  rate  and  defect 
density  in  terms  of  the  local  chemical  environment  at  the  substrate.  A  simple  two- 
parameter  expression  for  growth  rate  is  obtained,  which  with  suitable  parameter 
choices  reproduces  the  results  of  more  detailed  mechanisms  and  experiment  over  two 
orders  of  magnitude  in  growth  rate.  The  defect  formation  model  suggests  that  the 
achievable  growth  rate  at  specified  defect  density  scales  approximately  quadratically 
with  the  atomic  hydrogen  concentration  at  the  substrate. 


I.  INTRODUCTION 

The  chemical  vapor  deposition  (CVD)  of  polycrystalline  diamond  films  hats  ad¬ 
vanced  to  the  point  that  thick,  optically-transparent  films  with  thermal  conductivity 
near  that  of  natural  type  Ila  diamond  may  be  synthesized.  Due  to  its  unique  combina¬ 
tion  of  hardness,  optical  transparency,  high  thermal  conductivity,  and  high  bandgap, 
CVD  diamond  is  uniquely  suited  to  many  challenging  optical,  electronic,  and  ther¬ 
mal  management  applications.  Some  examples  currently  being  developed  include  heat 
spreaders  for  multi-chip  modules  [1],  cold  cathodes  [2],  infrared  optical  windows  [3], 
and  coatings  for  cutting  tools  [4].  However,  for  CVD  diamond  to  find  widespread  use 
in  these  applications,  costs  have  to  be  significantly  reduced  [5]. 

The  contributions  to  the  cost  of  CVD  diamond  include  both  fixed  costs  and  the 
costs  of  process  gases  and  energy.  To  minimize  the  fixed  costs  per  carat  of  diamond 
produced,  reactor  throughput  must  be  maximized,  which  requires  achieving  high  film 
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growth  rates.  To  minimize  the  gas  and  energy  costs,  the  proce.ss  efficiency  in  con¬ 
verting  raw  feedstocks  into  diamond  should  be  maximized.  Which  of  these  objectives 
is  most  important  will  depend  on  the  growth  process  employed  (for  example,  hot- 
filament,  microwave  plasma,  combustion  flame,  or  arcjet). 

In  either  case,  if  the  growth  conditions  are  altered  to  improve  growth  rate  or 
efficiency,  it  is  necessary  to  know  how  this  will  affect  film  quality.  Although  current 
theoretical  models  can  do  a  reasonable  job  of  predicting  growth  rates  for  specified 
process  conditions  [6,  7],  models  to  predict  quality  are  rare  [8- 1C]  and  are  largely 
untested. 

In  this  paper,  a  simplified  model  of  the  gas-surface  chemistry  occurring  during 
diamond  growth  is  developed.  Although  a  detailed  chemical  mechanism  of  diamond 
growth  would  be  complex,  involving  perhaps  hundreds  of  reactions,  the  essence  of 
diamond  growth  chemistry  can  be  reduced  to  a  few  generic  steps  (essentially  surface 
activation,  adsorption  on  open  sites,  desorption,  and  incorporation  into  the  lattice). 
By  considering  these  steps  alone,  we  arrive  at  a  simple  two-parameter  rate  law  for 
diamond  growth,  which  fits  mccisured  growth  rates  well  with  appropriate  parameter 
choices.  We  also  propose  a  simple  defect  formation  model,  which  leads  to  a  relation 
between  growth  rate,  defect  density  and  atomic  hydrogen  concentration  at  the  sub¬ 
strate.  While  no  data  currently  exist  against  which  to  compare  this  relation,  it  agrees 
qualitatively  with  experience  in  film  growth,  and  is  simple  enough  that  it  is  amenable 
to  experimental  test. 

One  conclusion  of  this  analysis  is  that  achieving  high  growth  rates  while  main¬ 
taining  film  quality  requires  a  high  concentration  of  atomic  hydrogen  at  the  surface. 
This  is  made  difficult  by  the  rapid  recombination  of  H  on  the  diamond  surface.  In 
a  companion  paper  to  this  one  [11],  the  issue  of  atomic  hydrogen  transport  to  the 
substrate  is  considered  in  depth.  Taken  together,  the  simplified  surface  chemistry 
models  presented  here  and  the  results  for  H  transport  in  Ref.  [1 1]  result  in  scaling  re¬ 
lations  for  growth  rate  and  defect  density  in  terms  of  controllable  process  parameters, 
and  allow  identification  of  the  optimal  conditions  for  high-rate  growth  of  high-quality 
diamond. 

Only  “conventional”  CVD  diamond  growth  chemistry  is  considered  here.  The 
feedstock  gas  is  assumed  to  be  a  dilute  hydrocarbon  in  a  carrier  gas  of  Hj,  and  the 
substrate  temperature  is  assumed  to  be  near  1200  K.  The  effects  of  oxygen  or  halogens 
on  diamond  growth  are  not  considered,  largely  because  the  chemistry  occurring  in 
these  systems  is  not  yet  well  understood.  Also,  the  effect  of  substrate  temperature, 
though  large,  is  beyond  the  scope  of  this  work.  Finally,  attention  is  confined  primarily 
to  the  high-quality  limit:  it  is  assumed  that  the  defect  density  is  small  enough  that 
it  is  not  necessary  to  account  for  large-scale  non-diamond  carbon  formation. 

In  the  sections  below,  we  consider  first  the  interactions  of  atomic  hydogen  with 
the  surface,  resulting  both  in  surface  activation  and  catalytic  recombination  of  H  to 
Ha,  followed  by  formulation  of  the  reduced  growth  mechanism  and  defect  formation 
model.  Finally,  these  models  are  used  to  classify  current  growth  processes. 


II.  SURFACE  ACTIVATION 


Diamond  growth  chemistry  is  initiated  by  abstraction  of  surface-terminating  hy¬ 
drogens  by  gas  phcise  atomic  hydrogen: 

CdH  +  H^Cj  +  Hj.  (1) 

Here  CdH  represents  a  generic  hydrogenated  surface  site  and  the  site  with  the 
hydrogen  removed  (a  radical  site).  Radical  sites  may  also  be  filled  again  with  H  via 

-h  H  CdH.  (2) 


Due  to  the  large  negative  change  in  free  energy  for  reaction  (2),  the  reverse  process 
(thermal  desorption  of  atomic  hydrogen)  is  negligible.  At  steady  state,  the  rate  of 
creation  of  radical  sites  by  reaction  (1)  just  balances  their  rate  of  destruction  by 
reactions  (-1)  and  (2).  Solving  for  the  steady  .state  radical  frarti'^p  /*  yield‘d 


r  =  ^  k,x„ 

^  [C5]  +  [CdH]  {k,+k,)XH  +  k.,' 


(3) 


where  ki  is  the  rate  constant  for  reaction  i,  and  is  the  atomic  hydrogen  mole 
fraction  at  the  surface.  (It  is  assumed  here  that  the  background  gas  is  primarily  H2. 
and  thus  [H]/[H2]  =  ''h-)  For  Xf{  sufficiently  large,  Eq.  (3)  shows  that  /*  reaches  a 
limiting  value  given  by 


/*  = 


fei 

+  ^2 ' 


(4) 


which  is  a  function  only  of  temperature.  Since  kijk^i  is  equal  to  the  equilibrium 
constant  A'p.i  for  reaction  (1),  equation  (4)  applies  as  long  as 


(5) 


To  estimate  fci  and  A:2,  most  studies  have  taken  the  approach  introduced  by  Fren- 
klach  and  Spear  [12],  in  which  diamond  surface  chemistry  is  assumed  to  be  similar 
to  the  chemistry  of  analogous  gas-phase  molecules  (alkanes).  A  typical  value  for 
abstraction  of  a  tertiary  hydrogen  from  an  alkane  is  [13] 

ki  «  1.3  X  10‘‘*  exp(— 7.3/Rr)  cm^/mol-s  (6) 


where  R  —  1.99  x  10~^  kcal/mol/K. 

The  recombination  rate  constant  ^2  niay  be  estimated  from  the  rate  constants  for 
H  recombination  with  alkyl  radicals.  Since  the  diamond  lattice  provides  a  means  to 
remove  the  enthalpy  of  recombination,  the  high-pressure  limit  of  the  gtis-phase  rate 
constants  is  most  relevant.  Data  on  H  recombination  with  alkyl  radicals  are  sparse. 
For  the  reaction  CH3  -{-  H  — ►  CH4,  Cobos  and  Troe  [14]  have  recently  concluded  that 
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the  high-pressure  recombination  rate  constant  is  2  x  10''*  cm^/mol-s,  independent  of 
temperature,  with  a  factor  of  2  uncertainty.  (Although  not  an  alkyl,  the  same  value 
is  reported  by  Ackermann  et  al.  [15]  for  H  recombination  with  the  benzyl  radical 
to  form  toluene.)  However,  for  H  recombination  with  the  isopropyl  radical  i-CsHr, 
Warnatz  [16]  recommends  a  high-pressure  rate  constant  of  2  x  10'^  cni'^/mol-s  -  an 
order  of  magnitude  lower.  Since  in  the  gas  phase  the  H  atom  may  react  with  either 
side  of  the  radical,  while  on  a  surface  only  one  side  is  accessible,  the  gas-phase  rate 
constant  should  be  divided  by  2  to  estimate  k2.  Therefore,  the  best  estimates  from 
gas-phase  data  are  that  ^2  is  in  the  range  10*^  -  10'*  cm^/mol-s. 

It  is  often  more  convenient  to  work  with  gas-surface  rate  constants  in  units  of  reac¬ 
tion  probability  per  site  per  collision,  rather  than  the  maiss-action  units  of  cm^/mol-s. 
(We  will  denote  the  rate  constant  in  reaction  probability  units  as  7;.)  Taking  reaction 
(1)  as  an  example,  the  conversion  from  k^  to  7,  may  be  done  as  follows.  The  proba¬ 
bility  per  unit  time  that  the  hydrogen  atom  on  a  given  CaH  site  will  be  abstracted  is 
simply  fci[H].  An  alternate  expression  for  the  abstraction  probability  per  unit  time  is 
where  71  is  the  abstraction  probability  per  site  per  collision,  Fw  is  the  flux 
of  atomic  hydrogen  striking  the  surface  (moles/cm^/s),  and  n,  is  the  total  surface  site 
density  (moles/cm^).  Since  F^  =  [H]c;^/4,  where  ch  is  the  mean  thermal  speed  of  an 
H  atom,  equating  these  two  expressions  for  the  abstraction  probability  leads  to 


7i 


4fcinj 

CH 


(7) 


Using  the  value  of  n,  appropriate  for  the  (Ill)  surface  (3  x  10  ^  moles/cm^),  Eq.  (6) 
is  equivalent  to 


7i 


100 

v/f 


exp(~7.3/«r). 


(8) 


At  1200  K,  this  expression  yields  71  ~  0.14. 

From  Eq.  (7),  at  1200  K  a  reaction  probability  per  collision  of  1.0  corresponds 
to  rate  constant  of  4  x  10'^  cm^/mol-s.  The  geis-phase  data  for  the  recombination 
reaction  would  then  suggest  that  72  should  be  at  least  0.25,  and  perhaps  >  1.0. 
(Note  that  values  greater  than  1.0  are  possible  at  low  radical  site  coverage,  since 
the  domain  of  attraction  for  the  reaction  may  exceed  the  distance  to  a  neighboring 
hydrogenated  site.)  It  is  possible,  however,  that  steric  hindrances  on  the  surface  may 
impede  recombination,  resulting  in  a  lower  value  for  72  than  would  be  expected  bcised 
on  the  gas-phase  data. 

An  alternate  approach  to  estimating  71  and  72  has  been  pursued  by  Brenner  and 
co-workers,  who  have  carried  out  molecular  dynamics  simulations  of  H  interacting 
with  the  (111)  surface  [17],  using  an  empirical  carbon-hydrogen  potential  [18].  This 
method  has  the  advantage  that  realistic  steric  factors  can  be  accounted  for,  eis  well  as 
energy  transfer  to  the  lattice.  The  value  for  71  they  obtain  (0.04  at  1200  K)  is  about 
a  factor  of  3  less  than  the  value  predicted  by  Eq.  (8).  For  72,  they  predict  a  value  of 
0.43  at  1200  K. 
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To  estimate  the  value  of  Xh  above  which  /’  reaches  its  limiting  value,  let  us  take 
^2/^1  ~  10?  based  on  the  results  of  Brenner  [17j.  The  best  estimate  to  date  of  A6’^2oo 
for  reaction  (1)  is  —16  kcal/mole,  calculated  for  the  ( 100)(2xl):lH  surface  [7],  (For 
other  surfaces,  AG®  should  be  within  about  2  kcal/mole  of  this  value.)  Therefore, 
A'p,i  ssj  800  at  1200  K.  Using  this  value  in  Eq.  (5),  /*  should  reach  the  saturation  value 
for  Xfi  above  about  1  x  10"'*.  Under  filament-assisted  growth  conditions,  the  atomic 
hydrogen  mole  fraction  has  been  measured  to  be  ~  1  x  10"^  at  the  surface  [19-22]. 
It  is  therefore  likely  that  /*  is  in  the  saturation  regime  under  these  conditions,  and 
increasing  the  atomic  hydrogen  mole  fraction  at  the  surface  will  not  result  in  greater 
surface  activation.  In  the  following  sections,  we  will  cissume  that  Xu  is  large  enough 
at  the  surface  that  /*  is  independent  of  gas  composition. 


III.  HETEROGENEOUS  RECOMBINATION  OF  ATOMIC  HYDROGEN 


In  addition  to  activating  the  diamond  surface,  reactions  (1)  and  (2)  constitute 
a  mechanism  for  surface-catalyzed  recombination  of  H  to  H2,  since  their  sun<  is  the 
reaction  2H  —*  H2.  Because  of  surface  recombination,  H  must  be  replenished  by 
diffusion  or  convection  to  maintain  a  steady-state  H  concentration  at  the  surface. 

The  surface  recombination  rate  is  specified  by  the  recombination  coefficient  7//, 
which  is  defined  as  the  probability  that  an  H  atom  recombines  when  striking  the 
surface.  The  H  surface  recombination  rate  is  therefore 

Rh  -  7//  ~r~  moles/cm*/s.  (9) 

4 

We  assume  that  recombination  is  first-order  [23],  and  therefore  7//  is  independent  of 
[Hj.  Assuming  the  only  two  types  of  surface  sites  present  are  CdH  and  Cj,  then  the 
recombination  coefficient  7//  is  simply 


7h  =  (1  -  Dli  +  /*72- 

Substituting  from  Eq.  (4)  for  /*,  this  reduces  to 

IH  2  V7i  72/ 


(10) 


(11) 


which  shows  that  the  recombination  coefficient  is  simply  the  harmonic  mean  of  the 
abstraction  and  recombination  probabilities.  Note  that  if  72  7i,  then  'yn  ~  27j. 

independent  of  72. 

Recently,  direct  measurements  of  jh  on  CVD  diamond  have  been  reported 
[23-25].  The  reported  values  are  shown  in  Fig.  1,  along  with  theoretical  estimates  us¬ 
ing  Eq.  (11).  The  mea.surements  by  Harris  and  Weiner  at  1200  K  and  by  Krasnoperov 
et  al.  aX  1 1 19  K  are  in  reasonable  agreement  with  the  molecular  dynamics  results,  and 
indicate  that  the  value  of  7//  at  1200  K  is  slightly  greater  than  0.1.  The  estimate  of 
7//  based  on  gas-phase  rates  appears  to  be  somewhat  too  large;  this  may  indicate  that 
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Eq.  (8)  overestimates  the  abstraction  rate  constant.  In  the  high-temperature  region, 
Krasnoperov  et  al.  report  a  measured  activation  energy  of  6.0  kcal/mole  [25].  Near 
room  temperature  there  is  more  scatter  in  the  experimental  data,  which  is  apparently 
due  to  sample- to-sample  variations  [26]. 

IV.  A  REDUCED  MECHANISM  FOR  DIAMOND  GROWTH 

There  is  now  considerable  evidence  that  the  dominant  precursor  to  diamond  is 
the  methyl  radical  CH3  [27-34],  although  diamond  may  grow  from  other  species 
well,  notably  acetylene  [35,  33,  32,  36],  with  lower  reaction  rate.  Several  elementary 
mechanisms  have  been  proposed  for  diamond  growth  both  from  methyl  and  acetylene 
[7,8, 12,37-42],  which  consider  in  detail  how  the  growth  monomer  (CH3  or  C2H2)  may 
be  added  to  a  particular  site  type  on  a  diamond  surface  through  a  series  of  elementary 
reactions. 

These  studies  have  demonstrated  that  chemically-reasonable  means  exist  to  add 
the  growth  species  to  the  lattice,  with  growth  rates  similar  to  those  measured  [39, 
7].  However,  they  do  not  attempt  to  provide  a  complete  description  of  diamond 
growth.  For  example,  little  work  has  been  done  on  elenieiitary  mechanisms  of  defect 
formation,  or  on  the  development  of  the  complex  surface  morphologies  seen  in  STM 
and  AFM  images  [43-48].  Due  to  the  complexity  of  these  problems,  progress  in 
applying  detailed  kinetic  mechanisms  to  issues  such  as  defect  formation  is  likely  to 
be  slow. 

Alternatively,  several  groups  [9,  10,  49]  have  begun  formulating  reduced  mech¬ 
anisms,  which  do  not  attempt  to  describe  each  step  in  detail,  but  seek  to  capture 
the  correct  qualitative  behavior.  In  this  way,  entire  classes  of  elementary  reactions, 
occurring  at  multiple  sites,  are  subsumed  into  a  single,  approximate  overall  reaction. 
By  keeping  the  number  of  overall  reactions  to  the  minimum  necessary,  it  is  possible  to 
derive  simple  analytical  expressions  for  the  growth  rate  and  defect  density,  which  may 
be  compared  to  experiment.  In  this  section,  a  simple  reduced  mechanism  is  devel¬ 
oped  which  incorporates  the  most  important  features  of  more  detailed  mechanisms. 
In  formulating  this  model,  the  simplest  possible  assumptions  are  made  at  each  step, 
consistent  with  current  knowledge  of  diamond  growth  chemistry. 

The  key  processes  that  a  reduced  mechanism  must  account  for  are: 

(a)  Establishment  of  a  steady-state  surface  radical  site  coverage  via  reactions  (1)  and 

(2). 

(b)  Attachment  of  reactive  hydrocarbon  species  (unsaturated  molecules  or  radicals) 
to  the  surface  at  these  sites 

c;  +  C„H„i.CdA.  (12) 
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(c)  Removal  back  to  the  gas-phase  of  the  surface  adsorbates,  either  by  thermal  des¬ 
orption 

CaA C„H„„  (13) 

or  attack  by  atomic  hydrogen  (etching) 

CaA -f  <jH Ca  +  products.  (14) 


(d)  Incorporation  of  the  adsorbate  into  the  diamond  lattice  (with  abstraction  of 
adsorbate  hydrogens  by  H) 

CaA-t-(m-l)H-^CaH-h(m-l)H2.  (15) 

While  all  proposed  mechanisms  contain  steps  similar  to  (a),  (b),  and  (d),  not  all 
allow  for  removal  of  adsorbates  [step  (c)].  In  fact,  this  step  is  critical  for  faceted  film 
growth.  For  smooth  crystal  surfaces,  incorporation  must  occur  primarily  at  steps, 
which  then  move  across  the  crystal  face.  This  suggests  that  adsorbates  which  adsorb 
at  terrace  sites  cannot  incorporate  at  these  sites  —  they  must  either  migrate  across 
the  terrace  to  steps  where  they  are  incorporated,  or  must  return  to  the  ga.s  phase. 
Since  surface  migration  is  expected  to  be  negligible  (due  both  to  the  large  C-C  bond 
strength  and  the  termination  of  the  surface  with  H)  [50],  the  observation  of  faceted 
films  indicates  that  return  of  adsorbates  back  to  the  gas  phase,  either  through  thermal 
desorption  or  etching,  must  occur  during  growth. 

The  above  steps  do  not  explicitly  include  mechanisms  for  non-diamond  growth. 
It  is  assumed  that  conditions  are  such  that  high-quality  diamond  s  being  grown,  and 
defect  incorporation  is  a  rare  event.  A  simple  model  for  defect  formation  will  be  given 
in  the  next  section. 

It  is  clear  that  the  incorporation  reactions  must  involve  atomic  hydrogen,  since  the 
remaining  hydrogens  on  the  adsorbate  must  be  removed  by  abstraction.  In  general, 
this  will  not  occur  in  a  single,  elementary  step,  and  thus  the  rate  for  this  step  may 
have  a  complicated  dependence  on  [Hj.  However,  for  simplicity  and  in  the  absence 
of  information  to  the  contrary,  the  incorporation  step  is  assumed  to  be  first-order  in 
[Hj.  Similarly,  it  is  assumed  that  the  etching  step  is  first-order  in  [Hj. 

Since  incorporation  occurs  primarily  at  steps,  k;  should  be  proportional  to  the 
step  density  /,.  It  is  assumed  here,  for  simplicity,  that  the  step  density  is  constant. 

With  these  assumptions,  a  simple  kinetic  expression  may  be  written  for  the  rate 
of  change  of  the  surface  concentration  of  adsorbed  hydrocarbons  [CaA]: 

=  ^4CnHj[C5]  -  ^a[CdA]  -  (k,  -I-  A:,)[CdA][H].  (16) 
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(17) 


At  steady  state,  d[CdA|/df  =  0,  and  thus 

MC^m]{C3l 


[CdA]  = 


The  incorporation  rate  Rc  in  moles/cm^/s  is  simply 

Rc  =  A:.(CdAj(H]. 


(18) 


Dividing  Rc  by  the  molar  density  of  diamond  nj  (0.2939  moles/cm^)  gives  the  linear 
growth  rate: 


G  = 


MCdA](H] 


nd 


Substituting  from  Eq.  (17)  for  [CdA]  results  in 


where 


and 


9i  = 


gir[C„H„,][H] 

92  +  [H] 

kikaTls 


{kc  +  ki)nd 


(19) 


(20) 


(21) 


92  - 


kg  +  ki 


(22) 


Equation  (20)  shows  that  there  exists  a  critical  H  concentration,  equal  to  g2, 
abov-  ’vhich  the  growth  rate  is  independent  of  [H]  (assuming  /*  is  independent  of 
Xu,  as  discussed  above).  The  reasons  for  this  behavior  can  be  seen  by  examining 
the  limiting  cases  [H]  <C  Qi  and  [H]  52-  Eor  [H]  -C  92,  the  steady  state  adsorbate 

concentration  is  determined  by  a  balance  between  adsorption  and  thermal  desorption, 
and  therefore  is  independent  of  [Hj.  Since  the  incorporation  rate  is  assumed  first- 
order  in  [H]  [Eq.  (18)],  the  growth  rate  is  proportional  to  [H]  in  this  limit.  In  the 
opposite  limit  ([H]  ^2),  incorporation  and  etching  are  fast  compared  to  thermal 

desorption,  such  that  the  steady  state  adsorbate  concentration  is  determined  by  the 
balance  between  adsorption  and  incorporation  and/or  etching.  In  this  limit,  [CjA]  is 
inversely  proportional  to  [H],  and  therefore  the  growth  rate  becomes  independent  of 
[H], 

In  previous  numerical  modeling  work  [6,7],  we  have  shown  that  a  methyl  growth 
mechanism  proposed  by  Harris  [39]  lead  to  predicted  growth  rates  within  a  factor 
of  2  both  for  low-rate  methods  (hot-filament)  and  high-rate  methods  (oxyacetylene 
torch  and  DC  arcjet).  While  the  good  absolute  agreement  with  experiment  is  fortu¬ 
itous,  the  reproduction  of  the  correct  scaling  of  the  growth  rate  with  gas  composition 
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is  more  significant.  This  mechanism  considers  reactions  which  convert  the  molecule 
bicyclononane  to  adamantane,  and  postulates  that  diamond  growth  on  the  (100)  sur¬ 
face,  where  sites  similar  to  the  opposing-hydrogen  site  on  bicyclononane  are  present, 
proceeds  in  a  similar  manner.  The  original  proposal  considered  the  ( 100)( lx  1  ):2I1 
surface,  but  more  recently  Harris  and  Goodwin  (7]  have  shown  that  this  same  mech¬ 
anism  is  the  rate-limiting  portion  of  a  more  complex  mechanism  for  growth  on  the 
(100)(2xl  ):1H  surface,  which  STM  images  (43)  show  is  the  predominant  reconstruc¬ 
tion  present  on  the  (100)  face  during  CV'^D  growth. 

Although  the  Harris  mechanism  consists  of  24  reactions  (actually  12  reversible 
reactions),  the  growth  rate  predicted  by  this  mechanism  can  be  fit  well  by  the  func¬ 
tional  form  of  Eq.  (20),  suggesting  that  the  mechanism  can  be  represented  reasona!)ly 
accurately  by  a  smaller  reduced  mechanism.  The  growth  rate  in  ^m/hr  predicted  by 
the  Harris  mechanism  at  1200  K  is  well-fit  by  the  relation 


G  =  9  X  10’7* 


[CH3][H] 

5  X  10-9  -f  [H]’ 


(23) 


where  the  concentrations  are  in  moles/cm^.  This  simple  expression  reproduces 
the  growth  rate  predicted  by  the  Har’-is  mechanism  to  within  10%  over  the  range 
3  X  10~*°moles/cm^  <  [H]  <  10“^  cc,  10""moles/cm^  <  [CH3]  <  10“'’moles/cm^, 
over  which  range  the  growth  rate  varies  from  less  than  0.1  //m/hr  to  more  than  7000 
/tm/hr.  Using  the  value  of  /*  predicted  by  this  mechanism  in  the  high-quality  limit 
when  the  adsorbed  methyl  coverage  is  small  (/"  =  0.2)  leads  to 


G  = 


1.8  X  10“ 


[CHaKH] 

5  X  10-9  +  [H]‘ 


(24) 


Whether  or  not  the  Harris  mechanism  provides  a  correct  detailed  description  of  di¬ 
amond  growth  kinetics,  equation  (24)  provides  a  simple  relation  which  empirically 
predicts  measured  growth  rates  with  reasonable  accuracy. 

For  use  with  computational  models,  it  is  convenient  to  have  a  simple  set  of  surface 
reactions  which  reproduces  the  growth  rate  of  Eq.  (24).  and  which  also  produces  the 
correct  recombination  coefficient  of  H  on  diamond.  Since  the  growth  rate  predicted  by 
the  reduced  mechanism  depends  only  on  the  parameters  gif‘  and  g2,  many  different 
sets  of  rate  constants  can  reproduce  Eq.  (24).  Therefore,  some  rate  constants  ( asi 
be  set  arbitrarily.  Shown  in  Table  I  is  one  particular  implementation  of  the  general 
reduced  mechanism  [reactions  (1),  (2),  (12)-(15)],  designed  to  reproduce  Eq.  (24) 
and  also  give  an  H  recombination  coefficient  7//  =  0.12.  For  reactions  involving  a 
gais-phase  reactant,  the  rate  constant  is  listed  both  in  miiss-action  units  and  in  units 
of  reaction  probability  per  collision  per  site. 

In  the  mechanism  shown  in  Table  I,  the  etching  rate  constant  is  arbitrarily  set 
to  zero,  and  kd  (reaction  s4)  is  chosen  to  give  an  adsorbate  lifetime  on  the  surface 
of  100  /IS  [7].  Also,  72  is  set  to  0.4,  to  match  molecular  dynamics  simulation  results 
[17j,  and  71  is  set  to  give  the  desired  H  recombination  coefficient.  The  rate  constants 


9 


k.  ki  (reactions  s3  and  s5,  respectively)  are  then  chosen  to  yield  gif‘  and  g2 
•  dues  appropriate  to  Eq.  (24).  Note  that  CaM  represents  an  adsorbed  methyl  group, 
and  Cd  represents  an  atom  of  bulk  diamond.  The  incorporation  step  is  represented 
by  the  two  elementary  reactions  (s5)  and  (s6).  This  is  done  to  maintain  the  proper 
stoichiometry  for  the  incorporation  step  (CaM  +  2H  — ♦  CaH  +  Cd  +  2H2)  while  keeping 
it  first-order  in  (Hj.  The  rate  for  (s6)  is  set  to  a  large  enough  value  that  (s.5)  is  always 
the  rate-limiting  step. 


V.  DEFECT  GENERATION 


The  challenge  in  scaling  CVD  diamond  processes  to  higher  rates  is  to  do  so 
while  maintaining  high  film  quality  (low  defect  density).  To  do  this  requires  v! 
understanding  of  how  defect  generation  depends  on  the  local  chemical  environment 
There  are  many  types  of  point  and  extended  defects  present  in  CVD  diamond  films 
(including  sp^  carbon,  substitutional  impurities,  vacancie.s,  interstitials,  dislocations, 
twin  planes,  etc.  [51]).  At  present,  there  is  little  quantitative  information  relating 
defect  densities  to  growth  conditions. 

Here  a  generic  model  of  defect  formation  is  developed,  which  results  in  a  simple 
scaling  law  for  film  quality.  This  model  will  clearly  not  apply  to  all  types  of  defects. 
Nevertheless,  its  predictions  are  in  qualitative  accord  with  experience  in  film  growth, 
and  may  be  relevant  in  particular  for  hydrogen  and  sp^  carbon  incorporation. 

The  basic  assumption  of  this  model  is  that  defects  are  generated  when  an  adsor¬ 
bate  reacts  with  a  nearby  adsorbate  before  it  is  fully  incorporated  into  the  lattice. 
For  example,  two  neighboring  adsorbed  methyl  groups  could  react  to  form  an  sp^ 
ethylene-like  group,  which  then  could  be  overgrown  to  lock  in  an  sp^  defect. 

The  rate  of  defect  generation  Riiff  (moles/cm^/s)  is  assumed  to  be  proportional 
to  the  number  of  adsorbate  pairs  on  the  surface.  Assuming  randomly  distributed 
adsorbates, 

fide/  =  ^de/[CdA]^.  (25) 

where  kj^/  is  the  (temperature-dependent)  rate  constant  for  defect  formation.  The 
defect  fraction  in  the  film  Xde/  is  given  by  the  defect  formation  rate  divided  by  the 
rate  of  sp^  carbon  incorporation: 


Xdef  = 


fi, 


de/ 


fi, 


C 


Substituting  for  [CaA]  from  (19)  and  using  fic  =  ridG,  results  in 


(26) 


Xdrf  = 


^dkdef\  G 

k;  jpip- 


At  constant  substrate  temperature,  then, 


(27) 
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{28) 


^def  ^ 


G 


Equation  (28)  provides  a  very  simple  expression  relating  film  quality,  growth  rat(\ 
and  the  atomic  hydrogen  concentration  at  the  substrate.  The  qualitative  predirtiuiis 
of  Eq.  (28)  —  that  film  quality  and  growth  rate  are  inversely  related,  and  that  in 
creasing  atomic  hydrogen  at  the  surface  improves  film  quality  —  are  in  accord  with 
experience  in  a  wide  variety  of  growth  environments.  However,  it  must  be  emphasized 
that  the  quantitative  validity  of  this  expression  is  unclear,  and  at  present  there  exist 
no  known  experimental  data  available  to  directly  test  Eq.  (28).  In  particular,  the 
exponent  in  the  denominator  depends  on  two  assumptions  —  that  defect  formation 
is  zero-order  in  [H]  [Eq.  (25)]  and  that  incorporation  into  the  lattice  is  first-order  in 
[H]  [Eq.  (18)].  A  less-specific  generalization  of  this  relation  would  be 


A  dfj  oc 


(29) 


where  n  must  be  determined  experimentally. 

Some  insight  into  the  appropriate  value  for  n  can  be  gained  by  examining  the 
predictions  of  Eq.  (29)  for  conventional  filament-assisted  diamond  growth.  In  this 
environment  at  20  Torr,  the  H  concentration  at  the  surface  is  approximately  5  x  10~'" 
moles/cm^  [19,20],  and  therefore  from  Eq.  (2-3)  the  growth  rate  is  proportional  to 
[CH3][H].  Equation  (29)  then  becomes 


XdeJ  OC 


ICH3I 


Under  these  conditions,  the  ga.s-phase  reaction 


(30) 


CH4  +  H  ^CH3  +  H2 


(31) 


is  rapid  and  near  partial  equilibrium  [20,52].  As  a  consequence.  [CH3]  is  approxi¬ 
mately  proportional  to  [CH4][H]/[H2],  which  implies  that 


XdeJ  OC 


[CH4] 

[Hz][H]"-2 


(32) 


for  typical  hot-filament  conditions.  Taking  the  value  n  -  2  leads  to  the  prediction 
that  the  defect  density  at  a  given  substrate  temperature  is  determined  primarily  by 
the  input  methane  fraction  —  a  result  in  accord  with  much  experience  with  filament - 
assisted  reactors.  Therefore,  while  the  theoretical  justification  for  Eq.  (28)  is  tenuous, 
it  at  least  leads  to  reasonable  predictions  for  filament-assisted  conditions. 
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VI.  A  PROCESS  MAP 


The  simplified  models  discussed  above  indicate  that  the  growth  rate  and  defect 
density,  for  a  given  temperature,  are  functions  only  of  the  H  and  CH3  concentrations 
at  the  substrate.  In  Fig.  2,  contours  of  constant  growth  rate,  calculated  from  Eq.  (24), 
are  shown,  along  with  contours  of  constant  G/(H]^  (normalized  to  hot-filament  con¬ 
ditions),  which  provides  a  measure  of  relative  defect  density.  Typical  concentration 
regimes  of  most  common  diamond  CV'D  methods  are  also  shown.  For  the  hot- filament 
and  microwave  plasma  cases,  the  concentrations  are  taken  from  the  measurements  of 
Hsu  [19-21].  For  the  other  processes  shown,  the  H  and  CH3  concentrations  are  esti¬ 
mates  based  on  numerical  simulations  [6,53]. 

Several  points  are  apparent  from  Fig.  2.  First  of  all,  it  is  seen  that  standard  low- 
pressure  (20  Torr)  hot-filament  and  microwave  plasma  systems  operate  in  a  regime  in 
which  the  growth  rate  is  dependent  on  [H],  as  discussed  above.  On  the  other  hand, 
arcjet  reactors  may  operate  well  into  the  .saturated  regime.  Here  the  growth  rate  is 
independent  of  [H],  but  film  quality  still  improves  with  increasing  [Hj.  .Atmospheric 
pressure  combustion  torches  and  RF  plasma  torches  operate  in  a  transitional  regime 
between  these  two  limits. 

The  results  of  Fig.  2  also  suggest  that  the  film  quality  is  similar  for  hot-filament, 
microwave,  and  combustion  torch  films  grown  at  an  acetylene  to  oxygen  ratio  /?  =  1.1 
(although  the  effect  of  oxygen  on  quality  is  not  considered  here).  In  general,  as  tiie 
CH3  concentration  is  increased  at  constant  [H],  both  growth  rate  and  defect  density 
increase.  These  results  indicate  that,  in  the  absence  of  oxygen  or  halogen  addition, 
the  highest  quality  diamond  should  be  grown  in  high-speed  plasma  jet  reactors,  due 
to  the  high  H  concentrations  at  the  surface  achievable  in  these  systems.  They  also 
suggest  that  it  should  be  possible  in  such  plasma  systems  to  use  gas  compositions 
significantly  ric'  r  in  hydrocarbons  than  is  possible  in  most  other  methods,  again  due 
to  the  large  [H]  at  the  surface. 

These  results  show  that  achieving  high  growth  rates  while  maintaining  film  quality 
requires  delivering  a  large  amount  of  atomic  hydrogen  to  the  substrate.  In  fact,  the 
H  concentration  at  the  surface  is  the  primary  factor  limiting  the  achievable  growth 
rate  for  a  specified  film  quality.  Once  [H]  is  fixed,  the  methyl  concentration  may  be 
easily  set  by  adjusting  the  methane  fraction  in  the  feed  gas  to  give  the  desired  growth 
rate  and/or  defect  density. 
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TABLES 


TABLE  I.  A  reduced  methyl  diamond  growth  mechanism. 


Reaction 

k 

7 

Notes 

(si)  CjH  +  H  — ►  Cj  +  H2 

2.9  X  10*2 

0.07 

b,e 

(s2)  C5  +  H-*CdH 

1.7  X  10*3 

0.4 

d,e 

(s3)  Q  +  CHa-CdM 

3.3  X  10*2 

0.31 

c,e 

(s4)  CdM-q  +  CHa 

1.0  X  10^ 

- 

a,e 

{s5)  CdM  +  H  ^  CdM*  +  Hj 

2.0  X  10*2 

0.05 

c,e 

(s6)  CdM*  +  H-CdH  +  Cd  +  H2 

fast 

fast 

Units  for  k:  cm^,  moles,  seconds. 

•Arbitrarily  set 

’^Chosen  to  give  7//  =  0.12 

‘^Chosen  to  reproduce  Eq.  (24). 

‘^Chosen  to  match  molecular  dynamics  simulation  results  [17] 
'Substrate  temperature  =  1200  K. 
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FIGURES 


FIG.  1.  Recombination  coefficient  of  atomic  hydrogen  on  the  diamond  surface.  Mea¬ 
surements:  Harris  and  Weiner  [23]  (♦);  Proudfit  and  Cappelli  [24]  (#);  Krasnoperov  et  al. 
[25]  (solid  line).  Theory:  estimate  for  the  (111)  surface  using  Eq.  (8)  for  71  and  assuming 
72  =  1.0  (dashed  line);  estimate  based  on  molecular  dynamics  simulation  results  for  71  and 
72  for  the  (111)  surface  [17]  (□). 


FIG.  2.  Contours  of  constant  growth  rate  from  Eq.  (24)  (solid  lines)  and  of  constant 
relative  defect  density  from  Eq.  (28)  (dashed  lines)  vs.  the  methyl  and  H  concentrations 
at  the  surface.  Typical  operating  regimes  for  some  common  diamond  CVD  processes  are 
also  shown  [19-21,6,53].  The  hot-filament  and  microwave  values  are  for  typicad  diffusion- 
dominated  reactors  at  20  Torr.  The  combustion  torch  vailues  are  for  an  acetylene-oxygen 
flame  at  1  atm  with  C2H2/O2  =  1.1.  The  pressures  for  the  RF  torch  and  DC  arcjet  are  1 
atm  and  220  Torr,  respectively. 
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CHg  concentration  [moles/cm^] 


